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ABSTRACT 

Type I X-ray bursts are very common in neutron star X-ray binaries, but no 
Type I burst has been seen in the dozen or so binaries in which the accreting 
compact star is too massive to be a neutron star and therefore is identified as 
a black hole candidate. We have carried out a global linear stability analysis of 
the accumulating fuel on the surface of a compact star to identify the conditions 
under which thermonuclear bursts are triggered. Our analysis, which improves 
on previous calculations, reproduces the gross observational trends of bursts in 
neutron star systems. It further shows that, if black hole candidates have surfaces, 
they would very likely exhibit instabilities similar to those that lead to Type I 
bursts on neutron stars. The lack of bursts in black hole candidates is thus 
significant, and indicates that these objects have event horizons. We discuss 
possible caveats to this conclusion. 

Subject headings: accretion — black hole physics — X-rays: binaries, bursts 



1. Introduction 

When gas accretes onto a neutron star (NS) in a low-mass X-ray binary (LMXB), it is 
compressed and heated as it accumulates on the surface, leading to thermonuclear reactions. 
In many NS LMXBs, the reactions occur unsteadily and cause Type I X-ray bursts (Grindlay 
et al. 1976). Type I bursts have been observed in a large number of NS LMXBs (see Lewin, 
van Paradijs & Taam 1993; Strohmayer, Swank & Zhang 1998 for reviews), and the theory 
of these bursts is relatively well understood (Hansen & van Horn 1975; Woosley & Taam 
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1976; Joss 1977; Taam & Picklum 1978; Paczynski 1982, hereafter P82; Fujimoto et al. 1984, 
1987; Fushiki & Lamb 1987, hereafter FL87; Taam, Woosley & Lamb 1996; Bildsten 1998). 

Among the over 100 LMXBs known in the Galaxy, roughly a dozen systems, most of 
them transient sources (Tanaka & Shibazaki 1996), have been identified as black hole (BH) 
candidates. In these BH LMXBs, dynamical measurements give mass estimates for the 
accreting stars greater than the likely maximum mass ~ 3Mq of a NS (Shapiro & Teukolsky 
1983, hereafter ST83; Narayan, Garcia & McChntock 2001; and references therein). 

No Type I burst has been seen in a BH LMXB, even though, as we show in this 
Letter, BH LMXBs ought to produce bursts as efficiently as NS LMXBs if the accreting 
BH candidates possess surfaces. The lack of bursts thus represents possible evidence for the 
presence of event horizons. In §2 of this paper, we examine the stability of nuclear burning 
on the surface of a compact star. In §3, we discuss the predictions of the model for accretion 
onto a IAMq NS and a lOM© BH candidate with a hypothetical surface. We conclude in 
§4 with a discussion. 



2. The Model 

We consider a compact spherical star of mass M and radius R, accreting gas steadily 
at a rate S (gcm^^s^^). In the local frame, the gravitational acceleration is g = GM{1 + 
z)/E?, where the redshift z is given hy I + z = {I - Rs/ R)-^''^, and Rs = 2GM/c^ is the 
Schwarzschild radius. We assume that the accreting material has mass fractions Xq, Yq and 
Zq — 1 — Xq — Yq, oi hydrogen, hehum and heavier elements (mostly CNO). 

We consider a layer of accreted material of surface density T^max sitting on top of a 
substrate of fully burnt material (X = Y = 0, Z = 1). (The particular composition of the 
substrate is unimportant since we apply the inner boundary condition at its top.) Since 
the physical thickness of the accreted layer is much less than the radius, we work in plane 
parallel geometry and take g to be independent of depth. We solve for the density p, the 
temperature T, the outgoing flux F, and the hydrogen, helium and heavy element fractions, 
X,Y,Z — 1 — X — Y,as functions of the column density E (0 < E < T,max)- 

The evolution equations for the gas in the layer are the standard equations of stellar 
physics. H- and He-burning give 

dX _ _eH_ dY^ _ eH_ _ eue 

where for our problem the total time derivative takes the form d/dt = d/dt + E9/9E. Here, 
^H,He s-re the respective nuclear energy generation rates, and E*jj jj^ are the corresponding 
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energy release per unit mass of H and He burned (P82). For eh-i we include the pp chain and 
the CNO cycle, including fast-CNO burning, saturated CNO burning, and electron capture 
reactions, as described in Mathews & Dietrich (1984) and Bildsten & Gumming (1998). 
Since we are not concerned with modelling the bursts themselves, and since our stability 
criterion does not depend on the detailed treatment of the deep crust, wc do not include 
proton captures onto heavier nuclei; Schatz et al. (1999) illustrate some of the consequences 
of the rp-process burning on accreting neutron stars. For He-burning, we include the triple-a 
reaction, but not pycnonuclear reactions (e.g., ST83). We do not correct the reaction rates to 
include screening (e.g. FL87) since we are concerned only with determining whether nuclear 
burning of H and He can proceed stably under given conditions; stable burning of H and He 
utilizes almost exclusively the reactions included. 

Hydrostatic equilibrium gives dP/dT. — g. For the pressure P we use the expressions 
given in P82 for the gas, radiation and degeneracy pressure, along with a correction when the 
degenerate electrons become relativistic. Radiative transfer gives dT/dT, = {3K,F)/{16aT^), 
where a is the Stefan-Boltzmann constant and wc write the opacity k in the form 1/k — 
^1 i^rad + ^l i^cond- Wc cmploy Ibcu's (1975) fitting functions for the radiative opacity Krad, and 
an analytical formula from Clayton (1968), suitably modified for relativistic electrons, for 
the conductive opacity Kcond'i the latter expression agrees well with more modern treatments 
(e.g. Heyl & Hernquist 2001). Finally, the energy equation gives 

r^ds , , OF ,. 

pT—^p{eH + €He)+Pg^: (2) 

where s is the entropy per unit mass. The above five equations form a closed set. 

We have four outer boundary conditions at the surface of the star (S = 0). Two of these 
are {X,Y) = (Xo,lo)- We obtain the third boundary condition by equating the accretion 
luminosity of the infalling gas. Lace = 47ri?^Sc^2;/(l + z), to blackbody emission from the 
surface: AwR^aT^^^ — Lace- This gives the surface temperature Tout- Then, using Tout and 
an assumed value of Fout, we solve for the surface density profile p(S) from the radiative 
transfer equation, thus obtaining the fourth boundary condition. 

At the base of the accreted layer we have an inner boundary condition. We assume that 
the temperature at the top of the substrate. Tin — T{T,max), is fixed. We examine several 
values of Tj„ for layers with T,rnax — 10^, 10^° and 10^^ g cm~^. Applying the boundary 
condition at T^max rather than deeper down is an approximation, but the error due to this 

is not serious. For the high surface densities we consider, the heat transfer is dominated by 
conduction (Heyl & Hernquist 2001), so the temperature gradient for S > S^ax is small. 
Moreover, we examine models for several values of Tin which further mitigates any error. 

The calculations proceed in two stages. First, we solve for the steady state profile of the 
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accretion layer by making the replacement d/dt tld/dE, d/dT, d/dE in the governing 
equations. This gives five ordinary differential equations with four outer boundary conditions 
and one inner boundary condition. The solution gives the profiles of the basic fiuid quantities: 
X(E), r(E), P(E)=^E, p(E), r(E). 

Having calculated the steady state structure of the accretion layer, we next check its 
stability. Various local stability criteria have been discussed in the literature (e.g., Bildsten 
1998), in which one considers the properties of the gas at a single depth. FL87 proposed a 
global criterion involving an integral over the entire layer. While an improvement, this ap- 
proach is still unsatisfactory since the authors assumed a constant temperature perturbation 
throughout the layer. 

We have carried out a full linear stability analysis of the accretion layer. We start 
with the steady state solution and assume that it is slightly perturbed, (^(E) Q{E) + 
(5'(E) exp(7t), where Q corresponds to each of our five variables, and the perturbations 
Q'(E) are taken to be small. We linearize the five equations described earher (three of which 
include time derivatives), apply the boundary conditions, and solve for the eigenvalue 7. We 
obtain a large number of solutions for 7 (technically, there is an infinite number since we are 
dealing with continuous functions), many of which have both a real and imaginary part. We 
consider the accretion layer to be unstable if any eigenvalue 7 has a real part (growth rate) 
greater than the characteristic accretion rate jacc — ^f^max- 

If the steady-state model is unstable according to the above criterion, it cannot accrete 
matter stably with the particular T,max and E. Whether this instability manifests itself as a 
Type I burst depends on how the burning flame, once ignited, envelopes the surface of the 
star. This is as yet an unsolved problem, although Spitkovsky, Levin & Ushomirsky (2002) 
present a possible solution. In the following we assume that the instability will grow rapidly 
to the nonlinear regime and that the system will indeed exhibit Type I bursts. 

3. Results 

Figure 1 shows results for solar composition material (Xq — 0.7, Yq — 0.27, Zq — 0.03) 
accreting on two kinds of compact stars. 
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Fig. 1. — Regions of instability, shown by dots, as a function of accretion luminosity and 
stellar radius. Top left: IAMq NS with a base temperature Tj„ = 10^-^ K. Top center: 
Tin = lO^K. Top right: T^^ = lO^-^ K. Bottom left: lOM© BH candidate with a surface, and 
a base temperature Tj„ = lO'''-^ K. Bottom center: Tj„ = 10''' K. Bottom right: Ti„ = 10^-^ K. 
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The top three panels correspond to an accreting NS of mass IAMq. We consider a 
range of accretion rates, parameterized by the ratio Lacc/LEdd, where we take the Eddington 
luminosity to be LEdd = 4:'kGMc/ Kes with k^s = 0.4 cm^g~^. We also consider a range of 
radii for the NS: log{R/ Rs) = 0.2 — 0.6 {R = 6.5 — 16 km). For each choice of Lacc/LEdd 
and R/Rs, we try three values of the surface mass density of the accreted layer: T^max = 
10^, 10^°, 10^^ gcm~^. If any of the three cases is unstable, i.e., if it has any eigenvalue 
with Re(7) > 7acc, then we claim that the system will exhibit instabilities that may grow 
into Type I bursts. 

The results presented correspond to three choices of the temperature at the base of the 
accreted layer: Tin = 10^'^, 10^ and 10' ^ K from left to right. In this Letter we are primarily 
interested in transient LMXBs. Since these sources have very low luminosities in quiescence 
{Lx < 10^^ ergs~^), the core temperatures of the NSs are likely to be Tj„ < 10^ K (Possenti 
et al. 2001; Brown, Bildsten & Chang 2002). 

The calculations shown in Fig. 1 indicate that NSs are unstable to bursts for a wide 
range of Tj„, but that the width of the instability strip (as a function of Lace) is less for higher 
values of Tj„. The reason for this is clear from the analysis of P82 who showed that when the 
flux escaping from the stellar core into the accretion layer increases (which happens when 
Tin increases), bursting behavior is restricted to a smaller range of E. 

In the case of transient NS LMXBs, we expect T^n to be < 10^ K. Fig. 1 indicates that 
these sources should be unstable to bursts for accretion luminosities up to Lace ~ O.SLEdd- 
The predicted luminosity limit is generally consistent with observations; the only NSs that 
are known not to burst are bright Z sources with Lace Lsdd (e.g., Matsuba ct al. 1995), 
and X-ray pulsars. Although the latter are significantly less luminous than Eddington, they 
accrete effectively at close to the Eddington rate since the accreting matter is channeled 
onto a small area on the NS surface by strong magnetic fields (see Lamb 2000 for a detailed 
discussion of this argument). 

Below Lace ~ O-SLEdd, the instability is initially of a mixed form in which a He-burning 
instability triggers a burst in which H and He both burn explosively. At lower luminosities, 
nearly all the H is burned steadily and the instability corresponds to a pure He burst. These 
results are consistent with previous work (e.g., Bildsten 1998). The calculations indicate 
that bursting behavior cuts off below an accretion luminosity Laec ~ lO^^'^L^dd. The cutoff 
is the result of the restriction T^rnax ^ 10^^ gcm~^ in our models. Systems with luminosities 
below our cutoff are still unstable to He bursts, but only at extremely high column density 
'^max- We have ignored such bursts since the recurrence time tree is too long to be of interest 
for observations of transient X-ray binaries. 
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For Tin = 10^-^ K, there is a second instability strip at low luminosities L^cc ~ 10 Lem- 
This strip corresponds to pure H bursts; the possibility of such bursts was first noted by 
Fujimoto et al. (1987). An interesting difference between the two instability strips is that 
the strip on the right generally has complex eigenvalues 7 for the unstable modes while that 
on the left has real eigenvalues. As explained above, the gap between the two strips is real 
if we restrict ourselves to reasonable values of T^max and tree- We are not aware of clear 
observational evidence for or against the gap, but such evidence could be searched for in 
future observations. 

The bottom three panels in Fig. 1 show results for a IOM0 BH candidate with a surface. 
The three panels correspond to different choices of the base temperature: Tj„ = 10^'^, 10^ and 
]^q6.5 from left to right. The particular choices of Tin are motivated by the extraordinarily 
low quiescent luminosities of many transient BH LMXBs {Lx < 10^^ ergs~^, Narayan et al. 
2001; a few BH LMXBs are brighter than this hmit, but even these are not hkely to have 
Tin ^ 10'^'^ K). We consider stellar radii in the range \og(R/ Rs) — 0.05 — 0.45, corresponding 
to R between (9/8)Rs and about 3Rs, a factor of nearly twenty in the surface gravity. The 
choice {9/8)Rs corresponds to the smallest radius in general relativity for an object whose 
density either decreases or remains level with increasing radius (ST83). 

The calculations indicate that BH candidates with surfaces arc not very different from 
NSs in their bursting behavior. Except for a modest rightward shift of the positions of the 
instability zones, the results in the lower three panels arc quite similar to those for a NS with 
a similarly low value of Tj„ (upper right panel). As in the case of NSs, bursts arc expected 
for BH candidates even for larger values of Tin than we have considered, except that the 
instability strip becomes narrower when Tin exceeds 10® K. We conclude that BH candidates 
are as prone to the instabihties that lead to Type I bursts as NSs are. The absence of bursts 
in BH LMXBs is thus highly significant and argues for the lack of surfaces in these systems. 

Although our analysis is for a spherical, non-rotating lOM© star, the results are not 
expected to differ greatly for a rotating object. The effective surface gravity and accretion 
rate vary by only a factor of 2-3 from pole to equator even for a maximally rotating oblate 
ellipsoid that is on the verge of the secular triaxial instability; furthermore, the dependence 
on compactness and equation of state is small (see Gondek-Rosinska & Gourgoulhon 2002; 
ST83). In comparison. Fig. 1 shows that bursts are present for a wide range of surface 
gravity g (factor of 20) and accretion rate E (factor of 100). The rather modest variation of 
g and S with latitude in a rotating star is thus not likely to have an effect. 

Because our model focuses only on the most important physical effects, and neglects 
many details, the exact positions of the instability strips in Fig. 1 may be uncertain at the 
level of say a factor of two in accretion luminosity. We believe, however, that the overall 
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pattern of instability we have computed is fairly robust. Also, the calculations we present 
here do not directly predict what kind of bursts are produced. We leave this for a more 
detailed paper, but briefly, we find that for NSs with Lace ~ ^-^LEddi bursts have durations 
of a few seconds and t^e^ ~ 10 hours. For lower Lace-, as is well-understood, bursts have 
longer tree and larger fluences. These results agree qualitatively with observations. 

4. Discussion 

It is clear from the theory of bursts (Bildsten 1998, and references therein) that burst- 
ing behavior is largely determined by the surface gravity g, the mass accretion rate, and 
the composition of the accreting material. Since these parameters arc similar in NS and 
BH systems, the bursting behavior of the two should be similar. Wc have quantified this 
argument with a global linear stability analysis which represents a technical advance over 
previous calculations. 

The results presented in Fig. 1 show that if BH candidates had surfaces they ought 
to experience thermonuclear instabilities as commonly as NSs do; by inference, they should 
have frequent Type 1 bursts. However, no BH candidate has exhibited Type I bursts. The 
most obvious explanation is that NSs have surfaces and BH candidates have event horizons 
(Menou 2001). If there is no surface, material cannot accumulate, and therefore cannot 
become hot or dense enough to trigger a thermonuclear explosion. 

Before we can claim that this "proves" the reality of the event horizon, more work is 
needed. First, we need to show that the model is able to reproduce the more detailed features 
of Type I bursts as observed in NS LMXBs. The statistics of burst durations and recurrence 
times (e.g., van Paradijs, Penninx & Lewin 1988) and the occurrence of oscillatory behavior 
in some systems (Rcvnivtscv ct al. 2000) ought to appear naturally in a realistic model. 
Also, any NS systems that burst when they should not by our calculations, or do not burst 
when they should according to the model, need to be explained. Second, the role of the inner 
boundary condition needs to be explored in detail. In Fig. 1, we see that different choices of 
the base temperature for a BH candidate give similar results. We have tried other boundary 
conditions, and also tried changing the composition of the accreting gas. In all cases we 
find that the accumulating layer is unstable to bursts for a wide range of luminosity, both 
in IAMq NSs and lOM© compact stars with surfaces. Finally, the difficult issue of flame 
propagation over the surface of the star once the instability has been triggered needs to be 
addressed (e.g.. Lamb 2000, Spitkovsky et al. 2002); the effect of rotation, for instance, is 
presently unclear. 
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We should caution that the analysis presented here assumes that the accumulating gas on 

the surface of a BH candidate behaves like normal matter, with nucleons and electrons; this is 
the case for standard models of strange stars and Q-stars (Glendenning 1997). Obviously, the 
argument becomes invalid if the properties of the gas change drastically, e.g., if the nuclei 
disappeared and were replaced by quarks. Whether such extreme changes are plausible 
remains to be seen. The density and pressure at the base of the bursting layer do not go 
above few x lO^gcm"^ and 10^^ ergcm"^ even in the most extreme cases we have considered. 
It is hard to imagine exotic physics being important under these conditions (Glendenning 
1997). 

On the observational front, we should check whether some NSs that lie within the 
unstable regions of Fig. 1 are stable to bursts. Any obvious large-scale disagreement between 
the observed burst behavior of NS LMXBs and the results presented here would indicate 
that the model is missing important physics. In the case of BH binaries, we should use 
observational data to derive quantitative limits on bursting activity. The transient BH 
LMXBs are particularly important for such work since they vary over a wide range of Lace 
during their accretion outbursts (Tanaka & Shibazaki 1996). If these sources have surfaces 
it is virtually impossible to arrange matters such that the objects have no bursts at all over 
the entire range of Lace- A firm demonstration that BH transients do not have Type I bursts 
would thus be a strong argument for the presence of event horizons in these systems. 
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